Improving relativistic MOND with Galileon k-mouflage 
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We propose a simple field theory reproducing the MOND phenomenology at galaxy scale, while 
predicting negligible deviations from general relativity at small scales thanks to an extended Vain- 
shtein ( "k-moufiage" ) mechanism induced by a covariant Galileon-type Lagrangian. The model 
passes solar-system tests at the post-Newtonian order, including those of local Lorentz invariance, 
and its anomalous forces in binary-pulsar systems are orders of magnitude smaller than the tightest 
experimental constraints. The large-distance behavior is obtained as in Bekenstein's tensor-vector- 
scalar (TeVeS) model, but with several simplifications. In particular, no fine-tuned function is 
needed to interpolate between the MOND and Newtonian regimes, and no dynamics needs to be 
defined for the vector field because preferred-frame effects are negligible at small distances. The field 
equations depend on second (and lower) derivatives, and avoid thus the generic instabilities related 
to higher derivatives. Their perturbative solution around a Schwarzschild background is remarkably 
simple to derive. We also underline why the proposed model is particularly efficient within the class 
of covariant Galileons. 
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Although the hypothesis of dark matter is consistent 
with a wide range of observations, it might be an ar- 
tifact of our Newtonian interpretation of experimental 
data if the gravitational 1 /r^ law happens not to be vahd 
at large distances. In 1983, Milgrom proposed a sim- 
ple phenomenological rule, caUed "modified Newtonian 
dynamics" (MOND) W\, depending on an acceleration 
scale uq: In the vicinity of a (baryonic) mass M , any test 
particle undergoes the standard Newtonian acceleration 
apf — GM/rc^ if ajq > ao, but a slower decreasing one 
a ~ y^aoojv — V GMuq jr if a < ao- This law happens to 
fit remarkably well galaxy rotation curves for a universal 
constant @] 



small distances, and not only its ratio to the gravitational 
one (as in the standard Vainshtcin mechanism). 

We consider a scalar-tensor theory of gravity defined 
by the action 



S = 
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where we use the sign convention of 13[ and notably the 
mostly plus signature. The physical metric g^J_y will be 
defined in Eq. ([7]) below. The scalar field if> is chosen 
dimensionless, and its kinetic term is the sum of the fol- 



(1) lowing three contributions, where s . 



Moreover, it automatically recovers the TuUy- Fisher 
law \^ and explains the correlation of dark matter pro- 
files with the baryonic ones Q- 

However, reproducing this phenomenological law in a 
consistent relativistic field theory happens to be quite dif- 
ficult, as illustrated by almost three decades of abundant 
literature. In the present paper, we shall refine on one of 
the best models proposed so far, in our opinion, namely 
the tensor-vector-scalar (TeVeS) theory constructed by 
Bekenstein and Sanders |5l4lOl| . We will show that a gen- 
eralized Galileon action allows us to suppress deviations 
from general relativity at small distances, thanks to an 
extension of the Vainshtein mechanism occurring in mas- 
sive gravity [ll|, ll2|. In the model we propose, the mag- 
nitude itself of the anomalous force tends towards zero at 



-^standard — 
•^MOND = 
■^Galileon — 
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Here e°'^^^ denotes the Levi-Civita tensor, related to the 
fully antisymmetrical symbol [a/375] (whose values are 
or ±1) by e"^T'5 = (_g,)-i/2 [ct^^J]. A small mass term 
— im^t/j^, with 1/m greater than the largest cluster sizes, 
might also be added to the above kinetic term. 

Denoting as tmond = yGM/oo the scale where 
MOND effects start to manifest around a galaxy of bary- 
onic mass M, and as ry = (8fcG'Mao)^^^/c < tmond 
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the Vainshtein radius below which scalar effects are sup- 
pressed, we will see below that -Cstandard dominates at 
very large distances r > tmond/^, ^mond at intermedi- 
ate ones ry < r < ryiom/^, and >CGaiiioon at small scales 
r < ry. As stated above, numerical fits of galaxy rotation 
curves 0] give the val ue ffl ) for ag, while their flatness 
up to about 10 rjviOND 14] imposes that the positive di- 
mensionless constant e is smaller than 0.1. The constant 
k entering /^GaUioon has dimension of a [length]'', and we 
will show below that a numerical value 
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(100 kpc)"* 



(6) 



allows the model to pass solar-system tests while pre- 
dicting MOND effects even for the lightest known dwarf 
galaxies. 

Several ingredients of action ([2]) are borrowed from the 
MOND literature. In particular, the scalar kinetic term 
-Cmond is well known to generate an extra acceleration 
y^GMoo/r on any test mass at distance r from a source 
of baryonic mass M 15| . Note that we did not write it 
simply in terms of |sp/^, but that an absolute value is 
involved only within the square root, to ensure that the 
scalar field carries positive energy whatever the sign of 
s [m . The standard kinetic term /^standard has a negligi- 
ble influence because it is multiplied by the small positive 
constant e, but it ensures that the dynamics of the scalar 
field is well defined when s passes through a vanishing 
value [ll,[l3- 

As in Refs. the matter action S'mattcr assumes 

that all matter fields ip are minimally coupled to a phys- 
ical metric 



= e ^'^g^^ - 2sinh(2(/7)C/^[/^ , 



(7) 



where is a timelike unit vector field, i.e., g^'^U^Ui, — 
— 1. Light deflection by galaxies and clusters would 
indeed be inconsistent with experiment if matter was 
merely coupled to the scalar fleld via a conformal metric 
.9^iy — e^'^g^iu [1, H, [ill- In a locally inertial frame where 
9iiu = diag(— 1, 1, 1, 1), and choosing the observer's ve- 
locity such that C/^ — (1,0,0,0) hes along his proper 
time direction, Eq. ([7]) merely means that goo = e^'^goo 
but Sij — e '^'^gij (mimicking thus the behavior of pure 
general relativity in presence of dark matter). Previ- 
ous fleld theories attempting at reproducing the MOND 
dynamics predicted large preferred-frame effects in the 
solar system, inconsistent with experiment, if was 
assumed to be a constant vector field 0, [lB|. This is 
whyit was assumed to be dynamical in the TeVeS model 
d, Q , with the idea that it could align with the mat- 
ter's local proper time direction. However, making the 
vector field dynamical by adding a kinetic term — i^^^, 
or any function of it, causes this vector to be unsta- 
ble [la, S, In the present paper, we will see that 
preferred-frame effects remain negligible in the solar sys- 
tem even if C/^ is assumed to be constant, thanks to the 



Vainshtein mechanism at small distances. This mecha- 
nism also allows us to choose other forms than ([7]), in 
contrast to Refs. We could also choose the physi- 

cal metric as « e^'^g^^ + B{ip, s)ip^^'p^v, but B would 
need to be a fine-tuned function of both Lp and its stan- 
dard kinetic term, s = (lyj^)^, in order to be consis- 
tent with the observed light deflection by galaxies and 
clusters. Moreover, the conditions for consistency of the 
theory within matter would be quite involved [16|]. The 
choice ([7|) , borrowed from [7l-[loj| , is thus the most natural 
one in the present framework. 

The reason why a mere action /^standard + -Cmond, 
Eqs. ([31) and ^ above, does not suffice to define a 
consistent relativistic field theory of MOND is that it 
would also predict an extra force ^jGMQao/r within 
the solar system (where Mq denotes the mass of the 
Sun), in addition to the Newtonian one GMq/t'^ and 
its post-Newtonian corrections. [We often use the word 
"force" instead of "acceleration" in the present paper, 
i.e., do not write the mass of the test particle to sim- 
plify.] This would be ruled out by tests of Kepler's 
third law and those of post-Newtonian dynamics. The 
literature considered thus "Relativistic AQUAdraric La- 
grangians" (RAQUAL), also known as "k-essence" theo- 
ries in the cosmological framework, i.e., a scalar kinetic 
term —\f{s) depending on a function of s = [ip^xj^. 
In order to reproduce the MOND dynamics, this kinetic 
term was assumed to take the form @ for small accel- 
erations (i.e., small values of s), while the scalar field 
behavior within the solar system depended on the shape 
of /(s) for large s. The clever choice of the literature 
was to impose that /'(s) tends towards a constant value 
for large s, say in order to recover a Brans-Dicke be- 
havior ip = —GMQ/f'^rc^ at small distances, so that the 
physical metric ([7]) reproduce the standard Schwarzschild 
solution up to a rescaling of the gravitational constant 
Gofi = G{1 + 1//^). The parametrized post-Newtonian 
(PPN) parameters /3 and 7 [131 take then strictly their 
general relativistic values /3 — j — 1, and classical solar- 
system tests are passed. However, binary-pulsar tests are 
directly sensitive to the matter-scalar coupling strength 
1 / , independently of the fact that the physical metric 
^ contains a "disformal" contribution proportional to 



UjjUui and they impose > 10* (2l|. As discussed in 
|ld |. such a large value is difficult to reconcile with the 
expression (|4]) needed for small s. One needs a fine-tuned 
interpolating function /(s), of a shape similar to Fig. 3 
of l6!|, in order to predict MOND effects while pass- 
ing binary-pulsar tests. Another idea would thus be to 
choose a function /(s) such that the scalar force at small 
distances is negligible with respect to the Newtonian 
one GMq jr^ , instead of keeping the same radial depen- 
dence ip'{r)(? = GMqI f^r^ . As underlined in [22], non- 
linear kinetic terms (i.e., those of k-essence/RAQUAL 
models) can reduce scalar effects at small distances, act- 
ing as a camouflage for the scalar (hence the name "k- 
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mouflage"). However, RAQUAL models must satisfy 
two conditions to have a Hamiltonian bounded by below 
and a well-posed Cauchy problem 16]: /'(s) > and 
2s/"(s) + /'(s) > 0. These conditions suffice to prove 
that ip"{r) < 0, i.e., that ip'{r) is a decreasing function of 
r, and the best we can obtain is thus an almost constant 
force ip' (r)c^ « oq within the solar system, the value ap 
being imposed by the MOND regime for r ~ tmond . But 
solar-system tests are precise enough to rule out a con- 
stant anomalous acceleration even numerically as tiny as 
([T]). We must therefore look for other possible scalar ki- 
netic terms, and this is where generalized Galileons enter 
our discussion. 

Galileons were first introduced in the cosmological con- 
text in |23| (a lthough they had actually already been 
studied in [24| in a quite different framework). In flat 
spacetime, they are theories whose field equations depend 
only on second derivatives of a scalar field, but not on 
their lower (Oth and 1st) nor on higher derivatives. One of 
their initial motivations was to generalize the key features 
of the decoupling limit of the DGP brane model p5| . 
which yields an equation of motion oc (Diys)^ — {^■,pi,vY 
for a scalar degree of freedom, pla^dn^ a crucial role 
in cosmology [26 



27|. References [2S 



( showed how 

to extend Galileon models to curved spacetime without 
introducing higher derivatives (while making now first 
derivatives also enter them). The same Lagrangians can 
be obtained in a suitable limit of brane models includ- 
ing Gauss-Bonnet-Lovelock densities [30| , or from dimen- 
sional reduction of such densities (sij. They can also be 
extended, in any dimension, to arbitrary p- forms possibly 
coupled to each other [s^] , or to general nonlinear models 
whose field equations depend on at most second deriva- 
tives [ss'l. Throughout this paper, we call Galileons this 
full class of models, although they go beyond the initial 
ones of 23 1 . 

The new ingredient of the present paper is the third 
scalar kinetic term of action ([2]), /^Gaiiioon, which strongly 
suppresses all scalar-field effects at small distances, as we 
will show below. This action (up to the factor — 8A:/3) 
was obtained in Ref. 3lJ for the first time by dimensional 
reduction of the Gauss-Bonnet-Lovelock density. It is 
written here in the compact form of l 29l| , and is also in the 
general classes considered in 32, |33| . It is easy to check 
that all field equations involve at most second derivatives, 
because of the antisymmetry of the Levi-Civita tensors 
entering ([5]). 

Although it is straightforward to write the full field 
equations deriving from action ([2]), for both ip and the 
metric, it will suffice in the present paper to obtain the 
perturbative solution for a static and spherically sym- 
metric in a Schwarzschild background ds^ ~ — (1 — 
rs/r)c^dt^ + dr'^/{l-rs/r) + r'^dVi^, where = 2GM/c^ 
denotes the Schwarzschild radius. It can be checked a 
posteriori that the backreaction of the scalar field on the 
metric has negligible effect as compared to present ex- 



perimental bounds. We will even assume that r ^ 
to simplify the expressions. More detail will be provided 
in a forthcoming publication [s^. Denoting as before 
if' = dip/dr, and integrating once the field equation for 
<y9, wc find 



2 / ^ ' « 

erp^-, 



(8) 



where the origin of the different terms is obvious from 
([3|)-([5]). The large number of field derivatives involved in 
'Ccaiiicon, Eq. ([l]), is rcspousiblc for the negative power 
of r in the first term of This is the central idea of 
the Vainshtein mechanism, since it makes this first term 
dominate at small distances. Imposing now that — )■ 
for 7' oo (otherwise p would diverge at infinity), we 
can immediately write the unique solution of which 
is a mere second-order polynomial equation for p': 




We thus easily recover the asymptotic Brans-Dicke be- 
havior p' K, GMjer^f? for very large distances r, the 
MOND regime p' ~ \JGM uq / rc^ at intermediate scales, 
but a small derivative p' ~ r/\/8k at small distances. 
Note that this small-distance behavior of p' does not 
depend at all on the mass M, which not only gener- 
ates the right-hand side (source) of Eq. dS]), but also 
the background Schwarzschild geometry entering the first 
term of ([8]) via the Riemann tensor of ([5]). This univer- 
sal small-distance behavior means that, paradoxically, 
any body generates strictly the same scalar force p'c^ 
in its vicinity. Figure [1] illustrates the three regimes of 
solution Neglecting its e contribution, the max- 

imum value (/?max = (GMao/2/c)^/*/2c is reached at 
ry — {SkGMaoY^^ /c, which defines thus a transition 
radius. 

In order to predict the MOND phenomenology in a 
galaxy of baryonic mass M, we need ry < fMOND, i-e.. 



k < GMc^i 



The lightest dwarf galaxies for which we 



have evidence for dark matter or MOND effects give thus 
an upper bound for the constant k. Draco- like dwarfs 
correspond to a baryonic mass M between 10^ and 10^ 
solar masses, but even tiny clusters of only 10'^ Mq seem 
to be dominated by dark matter [3^. Let us thus be 
conservative and impose k k, G{\Q^Mq)c^ /^icl^- This is 
the numerical value given in Eq. ([6|) above. 

We can now estimate the order of magnitude of scalar 
effects within the solar system. Since we chose a physi- 
cal metric of the disformal form ([7]), like in TeVeS, one 
could naively conclude that the PPN parameters /? and 
7 keep their general relativistic values f5 ~ j = 1. How- 
ever, there is now no meaning to define such parameters, 
because the PPN formalism assumes that no length scale 
enters the theory, and it needs the gravitational poten- 
tial to be oc 1/r at the Newtonian order. Therefore, the 
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PPN formalism 120] (while the terms corresponding to a2 
and vanish in the present model 0): we find that their 



FIG. 1: Behavior of the scalar force ip' (? in the three regimes 
of the theory: ip' oc r ai small distances (k-mouflage mecha- 
nism), ip' oc 1/r at intermediate range (MOND regime), and 
ip' oc 1 /r^ at large distances (asymptotic Brans-Dicke theory) . 



consequences of our anomalous potential (p oc r'^ + const, 
at small distances cannot be analyzed in the standard 
way. Let us thus merely compare the anomalous force 
ip'c^ Ri rc^/VSk it generates on a test particle, with the 
post-Newtonian forces ~ {GMQ/r^) x {rg/r) which have 
been precisely tested in the solar system. Since their ra- 
tio is (r^/rs)^/-\/2fc « (r/22AU)'' for the chosen value 
(lU of the constant k, the scalar effects are thus smaller 
than 10^^ times post-Newtonian ones at the Earth's or- 
bit, i.e., negligible with respect to the best experimental 
constraints. This ratio becomes even smaller for inner 
planets, and is < 10^^ for Mercury's orbit. On the other 
hand, scalar effects grow for outer planets, but they re- 
main 2 X 10~^ smaller than post-Newtonian forces at the 
orbit of Mars, and 3 x 10^^ at Jupiter's. This is consis- 
tent with the most precise planetary data. When con- 
sidering the Moon's orbit, the contribution of the Earth 
to the Riemann tensor entering ([S]) dominates over that 
of the Sun, and (p is almost spherically symmetric with 
respect to the Earth's center. Denoting now as r the 
distance to this center, we get again the universal be- 
havior (p'c^ ss r(? jsf^^ which must be compared to the 
Newtonian accelerations caused by both the Earth and 
the Sun and their post-Newtonian corrections. We find 
that scalar effects on the Moon's motion are 10~® smaller 
than post-Newtonian ones, i.e., four orders of magnitude 
smaller than the tightest experimental constraints de- 
rived from Lunar Laser Ranging [2^ . 

Preferred-frame effects can be estimated in a similar 
way. Assuming that the solar system is moving with a ve- 
locity w with respect to the preferred frame where f/^ — 
(1, 0, 0, 0), we compute the contributions proportional to 
in goo and to in goi , and their radial derivatives give 
us the magnitude of the anomalous scalar forces. Com- 
paring them to those generated by the ai-term in the 



ratio is < 8r3/(r'5ai V2fc) « (r/10"*AU)^/Q;i, giving thus 
scalar effects similar to those ofanai « 2x 10~^^ at the 
orbit of Mars. As above, when considering the Moon's or- 
bit around the Earth, the local value of pJ ~ r/\/8k must 
be used (r denoting now the Earth-Moon distance), and 
we find that preferred-frame effects caused by the scalar 
field are similar to those of an ai « 10^^^. Since the 
tightest constraint ai < 10~^ comes from Lunar Laser 
Ranging, we conclude that the model (H])-© does not 
predict any detectable violation of local Lorentz invari- 
ance in the solar system. 

Binary-pulsar tests are much more subtle to compute. 
A precise analysis would need either to study the (time- 
dependent) dynamics of a binary system at least up to 
order 0(I/c'^), or to be able to relate scalar multipoles at 
infinity to their local matter sources in spite of the non- 
linearities of the Vainshtein mechanism. In the present 
paper, we shall only estimate the rough order of mag- 
nitude of scalar radiation by comparing the local scalar 
forces to those of the precisely studied Brans-Dicke-like 
theories. Let us first note that the monopolar radia- 
tion, naively of order 0{l/c), is actually reduced to or- 
der 0{l/c^) because the local scalar solution generated 
by any body does not depend on time; this is similar 
to the case of standard scalar-tensor theories [s^l. On 
the other hand, the dipolar radiation starts at order 
0{l/c^) in spite of the fact that the local scalar solu- 
tion is the same around any body, independently of its 
mass. Indeed, if the two bodies of a binary system do 
not have the same mass (say, niA ^ rns), they do not 
move on the same orbit around their common center of 
mass, and the global scalar field they generate defines 
a preferred (oriented) direction in space. The dominant 
scalar radiation is thus a dipole, and we can estimate 
its order of magnitude by multiplying the one predicted 
in standard scalar-tensor theories (see, e.g. Eq. (6.52b) 
of Ref. [1^) by the square of the ratio of the present 
scalar force between the two bodies (rc^/-\/8fc) and the 
standard one. We get that the scalar field contribution 
to the time derivative of the orbital period is of order 
P ~ — GcP^{m A — mBYim A + rnB) I i^'^T^'^ km Am b), and 
numerically at least 10~^^ smaller than the tightest ex- 
perimental uncertainties. Although this estimate might 
be erroneous by some large numerical coefficient, we can 
anyway conclude that the present model should easily 
pass all binary-pulsar tests. 

Although the Galileon field equations involve at most 
second derivatives, and avoid thus the generic instability 
related to higher derivatives, this does not suffice to prove 
that these models are stable. One should carefully ana- 
lyze both the boundedness by below of their Hamiltonian 
density and the well-posedness of their Cauchy problem. 
The Hamiltonian of flat-space Galileons is straightfor- 
ward to derive 3J, |37[ , but the situation is much more 
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complex in curved spacetime, because all field equations 
for Lp and g^^ involve second derivatives of both of them. 
Around a given background, one should thus diagonalize 
the kinetic terms in order to test the stability of per- 
turbations and the hyperbolicity of their field equations. 
For instance, it would have no meaning to freeze g^j/ 
and study the perturbations of only in the scalar field 
equation (similarly to Brans-Dicke theory, which seems 
to contain a ghost scalar field for —3/2 < wbd < if one 
freezes the Jordan metric, whereas studying simultane- 
ously the dynamics of and Lp shows that the theory is 
stable even for such a slightly negative wbd)- This prob- 
lem of the consistency of Galileon field theories goes thus 
beyond the scope of this pap er, and we postpone it to a 
forthcoming publication |34| . 



Let us finally comment on our choice of Lagrangian ([5]) 
to obtain a k-mouflage mechanism reducing scalar effects 
at small distances. It happens to be the most efficient 
one amongst all those that we have analyzed. The 
highest-order Galileon Lagrangian which is nontrivial 
in 4-dimensional fiat space 23, 28, [2^, namely £5 = 

also generates a small (p'{r) w y/r / {^Ak^Y^'^ at small 
distances, but its slower radial dependence makes 
scalar effects still marginally detectable in the solar 
system (while giving again fully negligible preferred- 
frame effects). Imposing as above that the MOND 
phenomenology should occur in the lightest known 
dwarf galaxies, we get that scalar effects in the solar 
system are ~ (r/6AU)^/^ times post-Newtonian forces, 
i.e., about the order of magnitude of the most precise 
bounds. Similarly, scalar effects on the Moon's motion 
are 10~^ smaller than post- Newtonian ones, i.e., of 
the order of current limits. A full fit of planetary 
data taking into account the possible presence of such 
a small scalar force would thus be necessary to test 
whether it is already excluded or not. If not, this 
opens the exciting possibility to detect them in future 
higher-precision solar-system observations. All the 
other known covariant Galileon actions are either total 
derivatives in 4 dimensions, or yield a vanishing scalar 
field equation around a Schwarzschild background (like 
Eq. (23) of Ref. [3lj), or predict a negative p" {r) and 
therefore too large scalar forces in the solar system. 
On the other hand, generalized Galileon actions [33| 
involving non-differentiated fields and/or negative 
powers of s = (v?.a)^ can provide alternative models 34 1. 
but anyway less natural than ([5]). 



[1] M. Milgrom, Astrophys. J. 270, 365 (1983). 

[2] R. H. Sanders and S. S. McGaugh, Ann. Rev. Astron. 

Astrophys. 40, 263 (2002) arXiv:astro-ph/0204521 . 
[3] R. B. Tully and J. R. Fisher, Astron. Astrophys. 54, 661 



[4 

[6: 
[7: 

[9: 
[10: 

[11 

[12 
[13 
[14 

[15 

[16 

[17 

[18 
[19 

[20 

[21 

[22; 
[23; 
[24; 

[25; 

[26; 

[27 
[28; 
[29 
[30 
[31 
[32; 
[33 



Lett. 95, 171302 (2005) 
D 48, 3641 (1993) 
Astrophys. J. 429, 
480, 492 (1997) 



JHW2004, 012 (2005) 



(1977). 

S. S. McG augh, Phys. Rev 
[arXiv:astro -ph/0509305 . 
J. D. Bekenstein, Phys. Rev 
arXiv:gr-qc/9211017 . 
J. D. Bekenstein and R. H. Sanders 
480 (1994) ;arXiv:astro-ph/9311062 
R. H. Sanders, Astrophys. J. 
^arXiv:a stro-ph/9612099 . 

J. D. Bekenstein, Phys. Rev. D 70, 083509 

(2004) [Erratum-ibid. D 71, 069901 (2005)] 

arXiv:astro-ph/0403694 . 
J. D. Bekenstein, PoS 

arXiv:astro-ph/0412652 . 
R. H. Sanders, Mon. Not. Roy. Astron. Soc. 363, 459 

(2005) arXiv:astro-ph/0502222 . 

A. I. Vainshtein, Phys. Lett. B 39, 393 (1972). 

E. Babichev, C. Deffayet, and R. Ziour, Phys. Rev. D 82, 

104008 (2010) arXiv: 1007.4506 [gr-qc]]. 

C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravi- 
tation (Freeman, San Francisco, 1973). 

G. Gentile, P. Salucci, U. Klein, and G. L. Granato, 
Mon. Not. Roy Astron. Soc. 375 (2007) 199 
arXiv:astro-ph/0611355 . 

J. Bekenstein and M. Milgrom, Astrophys. J. 286, 7 
(1984). 

J.-P. Bruneton and G. Esposito-Farese, Phys. Rev. D 76, 

124012 (2007) arXiv:0705.4043 [gr-qc]]. 

E. Babichev, V. Mukhanov and A. Vikman, JHEP 0802, 

101 (2008) arXiv:0708.0561 [hep-th]]. 

M. A. Clayton (2001), arXiv:gr-qc/0104i03l 

G. Esposito-Farese, C. Pitrou, and J.-P. Uzan, Phys. Rev. 

D 81, 063519 (2010) ^arXiv:0912.0481 [gr-qc]]. 

G. M. Will, Living Rev. Rel. 9 3 (2006), 

< http: / /www.livingre views.org/lrr-2006-3 [ > , 

;arXiv:gr-qc/0510072 . 

G. Esposito-Farese, Proceedings of the 10th Marcel 
Grossmann Meeting (MGIO), World Scientific (2005) 647 
arXiv:gr-qc/0402007 . 

E. Babichev, G. Deffayet, and R. Ziour, Int. J. Mod. 
Phys. D 18, 2147 (2009) arXiv:0905.2943 [hep-th]]. 
A. Nicolis, R. Rattazzi, and E. Trincherini, Phys. Rev. D 
79, 064036 (2009) arXiv:0811.2197 [hep-th]]. 

D. B. Fairlie, J. Govaerts, and A. Morozov, Nucl. Phys. 
B373, 214 (1992) arXiv:hep-th/9110022^; D. B. Fair- 
lie and J. Govaerts, J. Math. Phys. 33, 3543 (1992) 
'arXiv:hep-th/9204074; . 

G. R. Dvali, G. Gabadadze, and M. Porrati, Phys. Lett. 

B 485, 208 (2000) arXiv:hep-th70005016 . 

C. Deffayet, G. R. Dvali and G. Gabadadze, Phys. Rev. 

D 65, 044023 (2002) arXiv:astro-ph /0105068 . 

C. Deffaye t, Phys. Lett. B 502, 199 (2001) 

[arXiv:hep- th/0010186 . 

C. Deffayet, G. Esposito-Farese, and A. Vikman, Phys. 
Rev. D 79, 084003 (2009) arXiv:0901.1314 [hep-th]]. 
C. Deffayet, S. Deser, and G. Esposito-Farese, Phys. Rev. 
D 80, 064015 (2009) arXiv:0906.1967 [gr-qc]]. 
C. de Rham and A. J. ToUey, JCAP 1005, 015 (2010) 
arXiv:1003.5917 [hep-th]]. 

K. Van Acoleyen and J. Van Doorsselaere, Phys. Rev. D 

83, 084025 (2011) arXiv:1102.0487 [gr-qc]]. 

C. Deffayet, S. Deser, and G. Esposito-Farese, Phys. Rev. 

D 82, 061501 (2010) arXiv: 1007.5278 [gr-qc]]. 

C. Deffayet, X. Gao, D. A. Steer, and G. Zahariade, 



6 



arXiv:1103.3260 [hep-th] . 
[34] E. Babichev, C. Deffayet, G. Esposito-Farese, 

D. A. Steer, and G. Zahariade, in preparation. 
[35] J. Wol f et al, Mon. Not. Roy. Astron. Soc. 406, 1220 

(2010) [arXiv:09Q 8.2995 [astro-ph.CO]]. 



[36] T. Damour and G. Esposito-Farese, Class. Quant. Grav. 
9, 2093 (1992). 

[37] S. Y. Zhou, Phys. Rev. D 83, 064005 (2011) 
[arXiv:1011.0863 [hep-th]]. 



